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Abstract 

A two-layer Heisenberg antiferromagnet is studied as a model of the bilayer 
cuprate YBa2Cu3C>6. Quantum Monte Carlo results are presented for the 
temperature dependence of the spin correlation length, the static structure 
factor, the magnetic susceptibility, and the 63 Cu NMR spin-echo decay rate 
1/^2G- As expected, when the ratio J2/J1 of the intrabilayer and in-plane 
coupling strengths is small, increasing J2 pushes the system deeper inside the 
renormalized classical regime. Even for J2/J1 as small as 0.1 the correlations 
are considerably enhanced at temperatures as high as T j J\ ~ 0.4 — 0.5. This 
has a significant effect on I/T^g, and it is suggested that measurements of 
this quantity at high temperatures can reveal the strength of the intrabilayer 
coupling in YBa2Cu306- 
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One of the unresolved issues regarding the high-T c cuprate superconductors is the role 
of the coupling between Cu02 planes within the same "block" in bi- and tri-layer com- 
pounds. Suggestions^ that this coupling might be responsible for the spin-gap behavior 
observed! in the bilayer cuprate YBa2Cu306.6 has spurred recent theoretical work on mod- 
els of two coupled Cu02 planes.! - ! Experimentally, the strength of the intrabilayer coupling 
in YBa2Cu306+x has not yet been accurately determined. Neutron scattering experiments 
performed in the insulating regime have not detected the high-energy branch of the spin- 
fluctuation spectrum up to energies of 60meV.@ This negative result analyzed in the frame- 
work of a linear spin-wave theory suggests a lower bound of J2 = 8 meV.i In the recently 
synthesized high-T c superconductor YBa 4 Cu70i 5 the two planes constituting a bilayer have 
slightly different chemical environments,! enabling various cross-relaxation experimentsE 
Analyzing such experiments by Stern et al.^ Millis and Monien estimated 5 meV < J2 
< 20 meV.i On the other hand, quantum chemical calculations indicate a value of J2 in 
YBa 2 Cu 3 6+ , as high as -50 meV.0 In this situation, it is important to consider new 
possible experiments that could accurately determine the magnitude of J 2 . 

Here we explore possibilities of extracting the strength of the intrabilayer coupling from 
experiments in the antiferromagnetic insulating regime, e.g. for YBa2Cu306- We report 
results of quantum Monte Carlo (QMC) simulations of a two-layer Heisenberg model defined 
by the hamiltonian 

H — J\ ^ ^ S a i ■ S a j + J2 X] ^1* ' ^2i, (1) 
tt=l,2 i 

where S a i is a spin-1/2 operator at site i of layer a, and denotes a pair of nearest- 

neighbor sites on a square lattice. This model should be a reasonable starting point for 
describing the magnetic properties of YBa2Cu30g at temperatures higher than the Neel 
temperature, which in clean samples is as high as 500K. Since the in-plane coupling J\ — 
1200-ff,! temperatures of interest are 0.4 < T / J\ < 0.5. We present results in this regime for 
the spin-spin correlation length, the static structure factor, the susceptibility, and the 63 Cu 
NMR spin-echo decay rate 1/T 2 g, and discuss the possibilities to experimentally detect the 



influence of a weak J 2- 

The finite-temperature physics of two-dimensional quantum antiferromagnets was ex- 
plained in the framework of a mapping onto a nonlinear a-model some time ago by 
Chakravarty, Halperin and Nelson.0 If the ground state is ordered the correlation length £ in 
the low-temperature renormalized classical (RC) regime diverges as e 27rps//T , where p s is the 
spin-stiffness constant. At temperatures T 3> p s , T is the only relevant energy scale and the 
behavior is quantum critical (QC), with £ ~ T _1 .0'E1 If p s is large, the cross over boundary 
between the RC and QC regimes of the cr-model may be at a temperature where a cr-model 
description of the antiferromagnet is no longer valid, and instead there is a direct cross-over 
from the RC to a local moment (LM) regime where the spins are essentially decoupled. 

For a single Heisenberg plane (J2/J1 = 0), the behavior at the temperatures of interest 
here is influenced by all the above regimes. For 0.4 < T / J\ ^ 0.6 there is a cross-over 
from RC to QC behavior, and at higher temperatures LM effects influence the behavior 
in the narrow QC regime J3il For the two-layer model described by ([I]), the ground state 
order increases with J 2 for small J2/J1, and has a maximum around J2/J1 ~ 0.8 before 
decreasing due to the tendency to singlet formation across the planes for larger J 2 . 10 A 
quantum phase transition to a disordered ground state occurs at {J2/ J\)c ~ 2. 5.00 The 
finite-temperature properties of near-critical systems have been studied numerically in detail, 
enabling a direct verification of the applicability of 1/N calculations for the nonlinear cr- 
model in the QC regime. Here we are concerned with values of J2I J\ more reasonable 
for modeling high-T c bilayer cuprates, and choose J2/J1 = 0.1,0.2, and 0.5. Owing to the 
enhanced ground state order for this range of J2/J1, one can expect the system to exhibit 
RC behavior at temperatures higher than for a single layer.! Below we present quantitative 
results for several experimentally accessible quantities. 

We begin by defining spin operators that are symmetric and antisymmetric with respect 
to interchange of the two layers: 

4« = s u ± S 2l . (2) 
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Using these we have calculated the correlation functions 



(3) 



the corresponding static structure factors 



(4) 



and static susceptibilities 




(5) 



where S±(i,r) =e rH S±(i)e~ rH , and L is the linear size of the system. The normalization 
in (f|) and ([|) has been chosen such that the standard definitions for a single plane are 
recovered for J 2 / 1 J\ = 0. For the numerical simulations we have used a QMC method based 
on stochastic series expansion!! (a generalization of Handscomb's methodill), which is free 
from errors of the "Trotter approximation" used in standard methods.il All results presented 
here are for 2x64x64 lattices (L = 64). 

In order to extract the correlation length we fit the correlation function C_(r) to the 
nonlinear cr-model forms discussed in Ref. [H5|. The results are shown in Fig. 1. For J 2 /Ji = 
we find good agreement with the results by Makivic and Ding.il The enhancement of the 
correlations with increasing J 2 /J± is evident. The correlation length quickly approaches the 
size of the lattice as T is lowered below T/Ji m 0.4 — 0.5, and in order to obtain reliable 
results for £ at lower temperatures lattices larger than L = 64 would have to be used. 

The correlation lengths graphed in Fig. 1 exhibit the exponential growth character- 
istic of the RC regime, as expected.il In order to more quantitatively address the ques- 
tion at what temperature the RC description is valid for a given J 2 /J± we study the ratio 
S-.(ir, 7r)/[Tx_(7r, 71*)]. In a classical system one alway has S/(xT) = 1 (for any q). A char- 
acteristic of the quantum antiferromagnet in the RC regime is that this relation remains 
satisfied close to the antiferromagnetic wave-vector. @ On the other hand, for a single plane 
in the QC regime S(ir, 7r)/[Tx(7r, 7r)] ~ l.lJU'illll For a two-layer model in the QC regime, 
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the gapped mode cannot be neglected if the temperature is of the same order as the gap, 
A oc v 'J1J2, and S^vr, 7r)/[Tx(7r, 7r)] will probably differ slightly from the above one-layer 
cr-model prediction. In any case, the RC value for this quantity should be 1, and is hence 
useful for determining whether the system is in the RC regime or not. The results displayed 
in Fig. 2 clearly indicate that the coupled planes approach the RC regime considerably faster 
than a single plane. 

The frequency integrated neutron scattering intensity is given by a combination of the 
structure factors S + (q) and S~(q), with weights depending on the momentum transfer per- 
pendicular to the planes.i Since the fluctuations corresponding to S + (q) are gapped at 
q = (tt, 7r), S + (ir,ir) saturates below a temperature set by the gap, whereas S-(n,7i) di- 
verges. Fig. 3 shows J 2 / J\ = 0.1 results for for q close to (ir, ir). The ratio S^(n, n)/S + (n, ir) 
is ~ 4 already at T/Ji = 0.5, and increases rapidly with decreasing T. At T '/ ' J\ = 0.4 the 
ratio is almost 30 (not shown in the figure). The smallness of S + (q) for q close to (7r,7r) 
may make the detection of the high-energy mode difficult, in particular because the even 
and odd modes cannot be completely separated experimentally.0 

In Fig. 4 shows results for the uniform susceptibility per spin, \ = x+(Q = 0)- The 
susceptibility for J 2 / J\ < 0.2 is very close to the single-plane result over the whole tempera- 
ture range considered here. The linear behavior seen for the single plane in the temperature 
regime 0.35 < T/J < 0.55 is in close quantitative agreement with the prediction for the QC 
regime,0 despite of the fact that the system actually crosses over to the RC regime at these 
temperatures (see Fig. 2).0@ An approximately linear behavior persists at these tempera- 
tures also for the systems with Jij J\ = 0.2 and 0.5, which are even deeper inside the RC 
regime. It would be interesting to compare the behavior with the QC and RC predictions for 
these couplings. As discussed above, this probably requires calculations for a full two-layer 
non-linear cx-model. 

We now turn to what we consider the most promising experiment for determining J2 
in the insulating regime. Recently, we showed^ that the NMR rates \jT\ and 1/T 2 g for 
La2Cu04 measured by Imai et a/S and Matsumura et ai.@ are well reproduced within the 



single-layer Heisenberg model and known hyperfine form factors (for the high-temperature 
regime, similar results were obtained by Sokol and co-workerslfE) . At high temperatures 
both 1/T\ and I/T2G show evidence of QC behavior, although the proximity to the RC 
and LM regimes influences the behavior as well. As shown above, even a small intrabilayer 
coupling pushes the system considerably further inside the RC regime. As a consequence, 
a system with J 2 /J\ ~ 0.1 — 0.2 should exhibit no QC behavior in the experimentally 
accessible temperature range. Below we present QMC results for I/T2G which should be 
useful for direct comparisons with experiments. 

The gaussian component of the spin-echo decay rate is related to the the nuclear spin- 
spin interactions mediated by the electronic spins. The coupling of a 63 Cu nuclear spin I 
at site in plane a to surrounding electronic spins is approximately given by the Mila-Rice 
formllEI 
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H = A ± (I*S: + I%S y a0 ) + A {l F S z a0 + BJ2Io- S a5 , (6) 

s 



where 5 denotes a nearest-neighbor of site 0. The hyperfine coupling constants A±, An and B 
are known from Knight shift measurements.il Pennington and Slichter derived the following 
form for l/T 2G , expected to be valid for the cupratesB 

1 rO.69 



ih 1 



E ^(M)l 1/2 ■ (7) 



T2G 

Here J z (xi,x 2 ) is the ^-component of the induced interaction between nuclei at X\ and x 2 



J z {x x , x 2 ) = -~ E M?\ - n)A(x 2 - rj)x(i - j), (8) 

where for the hyperfine coupling (^) one has A(0) = An, A(l) = B, and A{r) = otherwise. 
x(i — j) is the static susceptibility (H) written in coordinate space. Note that x in Eq. (|7|) 
runs over all the spins of both planes, except the spin at site in the plane where the nucleus 
considered resides. The constant 0.69 in ([?[) is the natural abundance of the 63 Cu isotope. 

In Fig. 5 we present results for J\/T 2 a versus T j J\ in units of K/s. The results expected 
in an experiment can be obtained by dividing with the relevant value of J\. We have 



used the standard experimental values for the hyperfine couplings; B = 41kOe//is and 
A\\ = -45.H In Ref. [L| we have shown that the QMC results for J 2 = agree well with the 
measurements on L^CuOj^S (the best agreement is obtained with a slightly smaller value 
for B; B « 37kOe//xs). We believe that our 1/T 2 g QMC results in Fig. 5 will be useful for 
determining the value of J 2 in YBa 2 Cu306, provided that measurements can be carried out 
in the regime of temperatures 500K < T < 800K. A potential difficulty is that J\ has to be 
known to rather high accuracy in order to establish the relation to the temperature scale of 
Fig. 5. 

In summary, we have studied the effects of a small intrabilayer coupling on the properties 
of the two-dimensional Heisenberg model. Even coupling ratios J 2 j J\ as small as 0.1-0.2 push 
the boundary of the RC regime up close to the highest temperatures accessible experimen- 
tally. This should have detectable consequences for a number of quantities. In particular, 
we suggest that measurements of the spin-echo decay rate 1/T 2 g at high temperatures would 
be useful for determining the strength of the intrabilayer coupling in YB^CuaOg. The spin- 
lattice relaxation rate, which we have not yet calculated for the two-layer model, should also 
be a sensitive probe. 
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FIGURES 



FIG. 1. The spin-spin correlation length vs. temperature for various strengths of the 
intrabilayer coupling. 

FIG. 2. The ratio S-(n, ir)/[Tx-(n, 7r)] vs. T for various intrabilayer couplings. In the 
RC regime this ratio is 1. The dashed line is the nonlinear cr-model result for a single plane. 

FIG. 3. The odd (upper panel) and even (lower panel) static structure factors for J 2 / ' J\ = 
0.1 close to the antiferromagnetic wave- vector q = (n, ir) at temperatures close to T/Ji = 0.5. 

FIG. 4. The temperature dependence of the uniform susceptibility per spin for various 
intrabilayer coupling strengths. 

Fig. 5. The 63 Cu spin-echo decay rate 1/T 2 g multiplied by the in-plane coupling J 1 vs. 
T/ Ji for various values of the intrabilayer coupling strength. 



10 




FIG. 1 




FIG. 2 




FIG. 3 



0.10 - 



0.08 - 



8 



O 
□ 



□ 



□ 



o □ 



0.06 - 



O J 2 /J 1= 0.0 
• J 2 /J 1= 0.1 
□ J 2 /J 1= 0.2 
■ ^=0.5 



0.04 



0.2 



0.4 0.6 

T/J, 



0.8 



FIG. 4 




FIG. 5 



